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ABSTRACT 



The general characteristics of buoyant momentum jets in water are 
described. Previous analytical modelling techniques utilizing the 
entrainment concept for prediction of trajectory and residual physical 
properties are discussed, and an overview of the existing experimental 
data base is given. The limitations of previous analytical modelling 
techniques are enumerated, generally resulting from incomplete or inade- 
quate equations of state. An existing comprehensive equation of state 
for pure and saline water is proposed for use in entrainment modelling. 
An original computerized procedure, based on appropriate conservation 
equations, is used to predict trajectory and physical properties of vari 
ous buoyant momentum jets. Comparison is made with previous analytical 
and experimental results for the cases of quiescent, flowing, and strati 
fied ambients. Finally, the comprehensive equation of state, coupled 
with the present computational procedure, is used to describe a com- 
plexly stratified ambient and the behavior of a buoyant momentum jet 
discharged into it. 
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NOMENCLATURE 



B - Characteristic jet width 

b - Dimensionless jet width, B/D 

c - Concentration 

D - Jet discharge diameter 

E - Volumetric entrainment 

p a ~p n 1/? 

F - Densimetric Froude number, UQ/(gD( )) ' 

2 p a" p m 

F, - Local densimetric Froude number, U / gB ( ) 

L "1 p o 

g - gravity 

Q - jet mass flow rate 

R - Ambient flow ratio, U /U n 

a u 

r - Radial jet coordinate 

S - Streamwise coordinate of jet veloci 

s - Dimensionless streamwise jet coordi 

t - Temperature 

U - Streamwise jet velocity 

U* - U ~ U cos 9, relative local veloci 

a 

u - Dimensionless streamwise jet veloci 

X - Horizontal Cartesian coordinate 

Z - Vertical Cartesian coordinate 

GREEK SYMBOLS 

a - Entrainment constant 

3 - Volumetric coefficient of thermal expansion 

y - Volumetric coefficient of concentration expansion 



ty 

nate, S/D; salinity 
ty 

ty, U/U Q 



o 



M ) - ( ), - ( )j 

9 - Local angle of inclination from horizontal 

X - Relative spreading ratio 
p - Density 
(j > - Azimuthal jet angle 

SUBSCRIPTS 
a - ambient 

e - at beginning of zone of established flow 
m - at jet centerline 
0 - at jet discharge 
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I. INTRODUCTION 



The cooling water discharge from a power plant into a large body of 
water, the thermally loaded condenser discharge from the condenser of a 
moving ship or submarine, and the high temperature gas issuing from a 
stack or gas turbine exhaust are all examples of buoyant momentum jets. 
The trajectory and behavior of such jets after discharge is influenced 
by factors such as initial jet velocity and buoyancy, ambient motion and 
stratification, and mixing rate. However, questions such as whether 
or not the jet will rise to a certain level, what the jet velocity and 
temperature will be at any point along its trajectory, or what effect 
ambient stratification will have on behavior, all require an involved 
quantitative analysis. Significant effort has been expended in the past 
few decades in attempting to understand the mechanics of buoyant jet 
mixing and trajectory, with the ultimate objective of developing accurate 
models to predict trajectory and decay. 

Certainly the need for such predictive models has grown. Contem- 
porary nuclear and fossil fueled power plants have thermal efficiencies 
on the order of 30-40%. The significant waste heat from these facilities 
takes the form of a thermally loaded discharge into either the atmosphere 
or a body of water. Sewage is often discharged as treated effluent into 
rivers, lakes, and oceans. The proper evaluation of the ecological 
impact of such discharges requires that their physical behavior be pre- 
dictable. More stringent environmental regulations and heightened public 
awareness place a premium on the accuracy of such prediction. 

The need to predict momentum jet behavior is not limited to environ- 
mental issues. Rapid advancement of the ability to detect small 
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temperature variations, concentration differentials, and turbulence 
anomalies may make it increasingly easy to detect various military craft 
and vehicles by virtue of propulsion system thermal discharges, wake 
turbulence, and wake concentration variations. The implications for 
weapons systems and platforms which rely on stealth for effectiveness 
are enormous. 

Given the wide range of applications in which an analysis of fluid jet 
behavior might be used, it becomes obvious that the range of possible jet 
or ambient characteristics that may be of interest is equally wide. 

Initial jet geometry, discharge parameters, degree of thermal loading, 
and turbulence characteristics, as well as ambient flow conditions, turbu- 
lence, and stratification, can be combined into an almost infinite number 
of scenarios. 

This investigation will be restricted to the case of a single, fully 
turbulent, circular water jet discharged into a surrounding water ambient. 
The case will be further restricted by limiting jet trajectory to two 
dimensions--in other words, ambient flow, if present, will be parallel 
to the horizontal component of jet velocity. Since jet encounter with 
an abrupt ambient discontinuity, such as a water-air interface, will not 
be addressed, the ambient will be considered to be infinite. 

Among the variables which will be considered are: 

(1) buoyancy effects, arising from density differentials between the 
jet and the surrounding ambient. These density differentials 
may arise from temperature and/or concentration variations. 

(2) ambient density stratification, arising from vertical non- 
uniformity of temperature and/or concentration in the ambient. 
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(3) ambient flow conditions, in which the magnitude and orientation 
of the ambient flow velocity relative to the jet are varied. 

(4) initial jet discharge characteristics, such as momentum 
orientation. 
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II. BASIC CHARACTERISTICS OF ROUND DISCHARGE SYSTEMS 



The terms "jet", "momentum jet", "forced plume" and "plume" are often 
used to qualitatively describe certain characteristics of a discharge 
system as it progresses through an ambient medium. It is generally under- 
stood that "jet", "momentum jet", and "forced plume" refer to that region 
where the momentum of the initial discharge is still sufficient to influ- 
ence jet behavior. "Plume" refers to a discharge in which the discharge 
momentum is either negligible to begin with, or small relative to the 
eventual total momentum produced by buoyancy. It is with this under- 
standing that these terms will be used. 

The jet/ambient system may be classified according to a number of 
characteristics: 

(1) Jet buoyancy 

(a) buoyant (positively or negatively) 

(b) neutrally buoyant 

(2) Initial jet orientation 

(a) horizontal (perpendicular to gravity field) 

(b) inclined 

(3) Ambient stratification 

(a) unstratified 

(b) linearly stratified 

(c) non-linearly stratified 

(4) Ambient motion 

(a) quiescent 

(b) flowing 
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Regardless of the classification of the jet or ambient, a jet passes 
through several flow regimes as it progresses along its trajectory. 

They are shown in Figure (2-1), and are: 

(1) The zone of flow establishment. In this region, flow character- 
istics are dominated by discharge conditions. Velocity and scalar 
quantity profiles (temperature, salinity, etc.) undergo transition 
from their initial discharge shapes through the action of a turbu- 
lent shear layer formed on the jet periphery. As mixing with the 
ambient progresses, the turbulent shear layer grows inward and 

the extent of the core of undisturbed profiles becomes smaller. The 
zone of flow establishment ends at the point where turbulent mixing 
reaches the jet centerline. The jet behavior in this region is 
strongly influenced by initial momentum and discharge conditions, 
and is only slightly influenced by the ambient. 

(2) The zone of established flow. This region begins when turbulent 
mixing reaches the jet centerline. The motion of the jet and its 
physical characteristics are governed by its momentum (initial and 
acquired), buoyancy, as well as ambient stratification and flow 
conditions. Initial discharge conditions play a progressively 
smaller role, and and the transition progresses from jet-like to 
plume-like behavior. 

(3) The far field. In this region the jet's initial momentum has 
negligible effect, and the jet may be convected by ambient flow. 

The jet fluid may be further diffused by ambient turbulence, and 

the distinction of the jet as a separate entity gradually disappears. 

The first two flow regimes constitute the near field, and will be the 
concern of this investigation. 
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Figure 2-1. Flow regimes and typical progression of a buoyant momentum jet 




Experimental work, begun by Albertson et al . [3] and continued and 
expanded by many others, has shown that within the zone of established 
flow, mean velocity profiles are nearly Gaussian: 

U = U m exp [-rW] (2-1) 

where is the mean centerline velocity, r is the radial jet coor- 
dinate, and B is a characteristic measure of jet width, or that radial 
distance at which U is equal to (1/e) of its mean centerline value. 

Profiles of jet scalar quantities, such as temperature and concen- 
tration, have also been found to be Gaussian in the zone of established 
flow by investigators such as Fan [7], Hoult et al . [16], and others. 
The profiles may be expressed as: 



At 



AC 



m 



t 



At exp|-r 2 /A 2 B 2 



f-rV 2 " 2 



Ac m exp |-r /A B 






( 2 - 2 ) 

(2-3) 



where At = (t-t ), At = (t-t,), Ac = (c-c ), Ac = (c -c ), and A is 
' a' m ' m a a' ’ m ' m a' 

the relative spreading ratio between velocity and density constituent 
scalar properties, a is defined as the inverse of the turbulent Schmidt 
Number. Figure (2-2) illustrates these profiles within the jet. 

A coordinate system to describe the trajectory and physical dimen- 
sions of a jet system is shown in Figure (2-3). The X coordinate is 
perpendicular to the gravity field and parallel to the flow (if any) of 
the ambient fluid. The Z coordinate is vertical and opposite to the 
gravity vector. The streamwise coordinate S defines the direction of 
mean centerline jet velocity at any point along its trajectory. The 
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Figure 2-2. Development of Gaussian velocity profiles 
in a momentum jet after discharge 
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Figure 2-3. Coordinate system for physical dimensions 
and trajectory of a jet system 
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local angle between S and X, or the inclination of the jet from the 
horizontal, is e. <p and r are polar coordinates defining the jet cross 
section, normal to S. Herein, any ambient medium motion is assumed to 
be horizontal. 

An important quantitative measure of relative momentum and buoyancy 
is the densimetric Froude number, F, given by 

U 0 

F = 2 (2-4) 

(gO ^0)1/2 

p 0 

The contribution of momentum is reflected in the numerator by the dis- 
charge velocity, Ug. The buoyancy effect is included in the denominator 
by the density differential term. Thus, the value of the densimetric 
Froude number ranges from near zero for plumes to infinity for pure, non' 
buoyant momentum jets. Hereafter, the term “Froude number" will be 
used to mean the densimetric form of Eqn. (2-4). 
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III. LITERATURE REVIEW 



A. REVIEW OF PAST MODELLING 

Several kinds of predictive models have been developed for the 
circular buoyant momentum jet. Although specific calculations consider 
different circumstances in origin of buoyancy, stratified/uniform ambi- 
ents, quiescent/coflowing ambients, etc., they all may be classified by 
basic method: 

(1) Algebraic models, based on either empirical data or simplification 
of differential model. These most typically predict only trajectory 
and jet width. Some, such as the model of Shirazi, McQuivey and 
Keefer [26] , also predict velocity, concentration, and temperature 
residuals. Data-based algebraic models tend to become unreliable 
when the basic conditions upon which they were based, such as 
general temperature and salinity range of the jet and ambient, are 
significantly changed. 

(2) Differential models, based on the relevant conservation equations 
(mass, momentum, energy, and scalar species). This modelling tech- 
nique allows prediction of jet trajectory and width, as well as 
velocity, temperature, and concentration decay downstream in the 
jet. Stratification and motion of the ambient may also be 
accommodated. 

Because of their limited scope, algebraic models will not be treated 
here. Certainly such models have a place in predictive use when the 
jet/ambient system involved is simple, and only information such as 
trajectory is required. However, the vast majority of effort in recent 
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years has involved the differential approach to jet modelling. In a 
majority of these differential models the entrainment mixing concept is 
invoked, rather than approaches utilizing mixing length hypotheses, k-e 
models, or eddy diffusivity. 

Morton et al . [20] were the first to use the entrainment concept to 
develop a buoyant jet model, as previously suggested by Taylor [27] . The 
concept supposes that the downstream induction of ambient fluid into the 
jet is proportional to the local jet centerline velocity, U m , and a 
characteristic jet width, B. Thus, 



E o< 2ttU B 
m 



where E represents volumetric rate of entrainment, or ambient inflow, 
into the jet, and is defined by 



dQ. . 

dS 



pE 



where Q is the total mass flow in the jet at any downstream location, s. 
Defining the constant of proportional ity , the entrainment constant or 
coefficient as a, the rate of entrainment can be written as: 



E = 2TTaU m B 

Solutions of the governing equations for differential modelling have 
been based on the following assumptions for round jets. 

(1) The jet flow is steady. 
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(2) The jet flow is fully turbulent. Molecular diffusion can be 
neglected in comparison with turbulent transport. 

(3) Streamwise turbulent transport is a negligible downstream trans- 
port mode, compared with streamwise convective transport. 

(4) Variation of fluid density throughout the flow field is small 
compared to a chosen reference density. Density variations are 
included only in buoyancy terms, the Boussinesq approximation. 

(5) Other fluid properties, such as viscosity, are constant over the 
range of interest. 

(6) Pressure is hydrostatic throughtout the flow field. 

(7) The jet remains axisymmetric throughout the near field. Velocity, 
temperature, density, and salinity profiles have no circumferen- 
tial dependence. 

The governing equations in the forms used in differential modelling 
are presented in Table (T3-1). 

With the exception of Hoult et al . [16], all studies cited in the 
following discussion have assumed velocity, temperature, salinity, and 
density profiles are Gaussian. This assumption, therefore, limits the 
applicability of such models to the zone of established flow. 

Hoult et al . circumvented the problem of having the model applicable 
only to the zone of established flow by assuming uniform ("top hat") 
profiles rather than a Gaussian distribution. This assumption was applied 
in the entire near field. As a result, the reduced form of the conserva- 
tion equations for this model differ from the reduced form used by others. 
Since Hoult's modelling technique is valid for both the zone of flow 
establishment as well as the zone of established flow, initial conditions 
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Table (T3-1). Governing Equations for Entrainment Modelling 
of Buoyant Momentum Jets 
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also differ. The end result of this approach is that values of p, t, c, 
and U ascribed to the jet at various points along the path are mean values 
for the entire jet cross section. This is a more limiting case than 
models using Gaussian profiles, where maximum values of jet properties 
result, and the entire cross section profile may be deduced from the 
appropriate Gaussian distribution. 

Abraham [2] initially used the vertical and horizontal momentum equa- 
tions, as well as the energy equation, to model jets discharged to a 
quiescent ambient. The continuity equation was not included. The 
solution required a pre-specification of the variation of B as a func- 
tion of S. Most other models have included the continuity equation in lieu 
of pre-speci fyi ng the B variation. 

The solution to the seven equations in Table (T3-1) yields values of 
jet centerline velocity, U , and temperature and concentration differ- 
ences At m and Ac m , as well as jet width D(s) and trajectory as functions 
of S. The solution of the equations, of course, also requires that the 
entrainment function E be specified. Herein lie the principal differ- 
ences between entrainment models. The models fall into two general cate- 
gories: those for a quiescent ambient, and those for a flowing ambient. 

1 . Quiescent Ambient Media 

Albertson et al . [3] and others have verified through measurements 
that for non-buoyant momentum jets, F = «, the appropriate value of 
within the zone of established flow is 0.057. There seems to be little 
disagreement with this value, based on numerous comparisons of differen- 
tial modelling and experimental data. 

Abraham [2] suggested, also on the basis of experimental evidence, 
that for relatively buoyant flows (small F), a = 0.085. This is in good 
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agreement with the suggestion of List and Imberger [18] of a = 0.082 
for pure buoyant plumes (F = 0). Fan [9] also suggested a = 0.082 for 
all flows except pure momentum jets. Fan also recommended, on the basis 
of his experiments, a = 0.057 for the pure momentum jet. 

In application, however, discharges are seldom either pure jets 
or plumes. Typically they are in some stage of transition away from 
jet behavior toward plume behavior. Morton et al . [20] proposed to model 
this transition by: 



a 2 

a = 0.057 + p— 

1 

where a 2 is an empirically determined coefficient, and F L is a local 
Froude number, based on the local density difference. The same general 
form was derived by Fox [10] for a vertically discharged buoyant jet. 

Hirst [14] maintained that for a discharge into a quiescent 
ambient, the entrainment function should depend on: 

(1) local mean flow conditions in the jet, i.e., U m and B. 

(2) local buoyancy within the jet, as indicated by local Froude 
number, and 

(3) jet orientation, 6 q. 

The following form was proposed: 

a = 0.057 + sin(e) 
r L 

This is the general form suggested by Morton and Fox, with the constants 
defined by fitting the function to known discharge and endpoint conditions 
of jet flow. 
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Another entrainment function for initially horizontal buoyant 
momentum jets is the jet-plume extrema fit proposed by Riester et al . 
[23]: 



a = jj0.057 cos e) 2 + (0.082 sin 0)^j 1/2 

From data on a buoyant jet discharged vertically downward into 
a quiescent ambient, Davis et al . [6] proposed: 



= 0.057 + 



0.083 

r .3 



A tabular summary of entrainment functions for discharges into quiescent 
ambient media is given in Table (T3-2). 

2. Flowing Ambient Media 

Hirst [15] proposed an entrainment function applicable to 3- 
dimensional buoyant jet flow, in which the horizontal component of initial 
jet velocity was not necessarily parallel to ambient flow. Eliminating 
such terms, the form for two dimensional buoyant jets becomes: 



a 2 

E = (a 1 + — sin 8)b[|U m -U cose| + a^sin 9 ] 



where the term 



! U-U cos 9 1 
1 m a 

represents the relative velocity of the jet with respect to the ambient, 
in the direction of jet flow. It is a pure "coflow" term. The term 
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Table (T3-2). Entrainment Functions for Use with Quiescent Ambients 
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(1) Applicable to simple momentum jet; (2) Applicable to simply buoyant plume; 

(3) Applicable to buoyant jet discharged at varying angles to a quiescent ambient; 

(4) Empirical fit for a buoyant jet discharged vertically downward to a quiescent ambient; 

(5) Applicable to buoyant jet discharged horizontally to a quiescent ambient; 

(6) Applicable to simple buoyant plume; ( 7)*Empi rical ly determined coefficient, vertically discharged. 



a, II sin 9 

O a 



represents the contribution to the total value of entrainment of pure 
cross flow, or ambient motion normal to the jet axis, as the jet turns 
due to buoyancy. 

Hirst specified values of a-j and so that the flowing ambient 
entrainment function reduced to the quiescent ambient function if II = 0. 

a 

He specified a. value of 9.0 for a^, based on a best fit to available data. 
This entrainment function thus became: 



E = (0.057 + %^sine)b[ |U-U cose| + 9.0U sin 6] 

r m a a 



Ginsberg and Ades [12] performed a least squares analysis on a 
large set of laboratory trajectory data. Using Hirst's entrainment model 
with a-j and a^ as specified, they found that a large variation in the 
value of a^ was necessary to fit predicted results with the data. They 
constructed a correlation for a^ as a function of F and the coflow ratio: 



a 3 = 25. 81 0[ F°* 19464 R 0,35155 ] . iq.825 



where the coflow ratio, R, is defined as: 



Other entrainment functions for coflowing ambients have appeared 
using variations of the governing differential equations. Schatzmann [25] 
proposed an entrainment function similar in form to that of Hirst, but 
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for use in a set of governing equations in which the Boussinesq approxi- 
mation was not invoked. Fan [9] proposed an entrainment function for 
coflowing ambients which included a drag term as well as the standard 
proportionality of entrainment with centerline velocity and jet width. 

He found that the value of the drag coefficient used, as well as the 
entrainment constant (a = 0.082 in this case) had to be readjusted to 
make the prediction conform with data with each change in discharge or 
ambient conditions. The entrainment functions of Hirst, Ginsberg and 
Ades, and Schatzmann are collected in Table (T3-3). 

The progression of differential modelling, once the entrainment 
function is specified, is similar in most cases. First, the governing 
equations are integrated over the jet cross section so that they appear 
in differential form. They are then non-dimensional ized with respect to 
chosen reference variables. Initial conditions are specified. The re- 
sulting equations are then numerically integrated over the desired range 
of the streamwise pathlength, S. However, different downstream trajec- 
tory and decay may be calculated for a given jet due to differences in 
any of the following: 

(1) Entrainment function chosen; 

(2) Initial conditions specified for the beginning of the zone of 
established flow; 

(3) Equation of state specified for the density of the fluid; 

(4) Computational technique. 

Examples of predicted jet trajectories and phsyical properties 
are presented in Figures (6-1) to (6-53). Comparison between models for 
identical or closely similar discharge and ambient conditions are pre- 
sented in Figures (6-11) to (6-15). Discussion of these predictions and 
comparisons will be made later. 
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With Coflowing Ambients: 
Complete Volumetric Entrainment 
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(8) Applicable to buoyant jet discharged at varying angles to a flowing ambient. Cross flow 
terms omitted. 

(9) Applicable to buoyant jet discharged at varying angles to a flowing ambient. 

(10) Same as [9], empirically determined coefficients 



B. REVIEW OF EXPERIMENTAL STUDIES 



By far the most comprehensive and often cited set of data for buoyant 
water jets resulted from the work of Fan [7]. The experiments concerned 
two classed of buoyant jets: 

(1) inclined jets discharged into a stagnant environment with linear 
density stratification. 

(2) buoyant jets discharged into a uniform cross stream, e Q = 90°, 
with no ambient stratification. 

The experiments for the flows of group (1) above were conducted in 
a 2.26m x 1.07m tank with a depth of 0.61m. The tank was stratified 
with successive 3 cm to 5 cm layers of aqueous salt solutions. Tank 
temperature remained constant within a 2°C range over the duration of the 
experiments. 

Nozzle diameters varied between 0.223 and 0.762 cm. Flow into the 
nozzle was provided from an unregulated head tank, which provided a 
discharge rate estimated by Fan to be constant within 3%. Measurements 
in each experimental run were limited to jet trajectory and half width, 
observed photographically by use of a tracer dye premixed into the dis- 
charged fluid. 

Fan described one of the dilemmas of conducting experiments on jets 
of small physical scale: 

For complete experimental check on theory, it is necessary to 
determine values of the jet velocity and density. Practically 
however, laboratory experiments on density stratified flows are 
usually limited in scale and do not allow the time required 
in measuring time fluctuating quantities. On the other hand, 
the jet trajectories and half widths can be determined con- 
veniently by photographic means. These two quantities are 
interrelated with other jet characteristics. Thus the com- 
parison of the observed and calculated values of these two 
quantities is believed to be indicative of the applicability 
of the theoretical solutions. [8] 
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The experimental discharge Froude number ranged from 10 to 60, with the 
exception of 3 runs with a non-buoyant momentum jet (F = °°) . 

In Fan's second group of experiments, with a flowing, unstratified 
pure water ambient, the saline jets were actually negatively buoyant. 
Conductivity measurements were taken by variable position probes at a 
number of downstream stations. The locus of stations of maximum concen- 
tration was defined as the jet centerline. The jet width was defined 
from concentration readings taken radially outward from the centerline. 

These experiments were conducted in a 40m flume, 1.1m wide, with a 
water depth of .51m. Flow was induced by inclining the flume. A region 
in the core of the flume flow with the least shear effects from wall 
boundaries was selected to introduce the jet via a nozzle. Variation of 
ambient flow in this region was estimated at +6% to -9%. 

Experimental runs were made for a Froude number range of 10 to 80, 
ambient flow of R = 0.0625 to 0.25, and a discharge diameter of 0.5 cm 
to 0.762 cm. 

In an attempt to ascertain the effect of ambient turbulence and shear 
introduced by the restricted cross-sectional dimensions of the flume, 
a limited set of concentration measurements were made in which the fluid 
in the flume was stagnant and the jet discharge towed through the ambient 
by a carriage mounted over the flume. In this case, conductivity probes 
were fixed to the carriage and moved to different relative positions in 
successive runs. 

The experiments of Fan have been used by numerous modellers, including 
Fan himself, as a basis for analytical/experimental comparison. The 
experiments also served to validate earlier hypotheses and observations 
regarding the Gaussian distribution of concentration in the zone of 
established flow. 
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Riester, Bajura, and Schwartz [23] studied horizontally discharged 
buoyant fresh and salt water jets. The ambient was quiescent and un- 
stratified. A 6.2 x 1.1 x 0.8 m tank was used in conjunction with a 0.87 
cm diameter discharge nozzle. Jet trajectory and width were recorded 
photographically by a tracer dye in the jet fluid. Temperature distri- 
butions were measured by a rake of thermocouples. Jet centerline was 
determined through the measured temperature distributions. 

A novel aspect of these experiments was the wide range of ambient 
temperatures util i ed (4.5°C-43.0°C) , and the use of both salt and fresh 
water jets. Some possible implications of the results will be discussed 
in a later section. 

Davis, Shirazi, and Siegel [6] measured the behavior of single and 
multiple port salt water discharges directed vertically downward into a 
fresh water ambient. A flowing ambient was simulated by mounting the 
discharge nozzle on a moving carriage mounted over a 17.1m tank. Water 
depth in the tank was 0.91m. Concentration profiles were measured by 
conductivity probes, and velocity measurements were made by hot film 
anemometry. Several runs of these experiments made with the carriage 
stationary led to the proposal of the entrainment function: 



= 0.057 + 



0.083 

-3 



These same experiments verified the Gaussian nature of velocity profiles 
in the zone of established flow. 

Shirazi, McQuivey, and Keefer [26] studied buoyant jets in flowing 
turbulent ambients. An inclined 120 ft. flume was used to produce 
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ambient flow. Turbulence was introduced by imposing a layer of varying 
sized rock in the flume bed. Temperature and salinity concentration 
were measured using conductivity probes. Turbulence was monitored by 
hot film anemometers. 

As a result of these experiments, a set of algebraic correlations 
were made, expressing temperature and concentration residuals, jet width, 
and trajectory as functions of downstream distance, F, and R. It was 
determined that as the level of turbulence increases, the rate of decay 
of centerline temperature and concentration with respect to the stream- 
wise coordinate increases. Not surprisi ngly , it was also found that 
correlation of the data became more difficult as turbulence level 
increased. 

Pryputniewicz and Bowley [22] conducted measurements of a buoyant 
jet discharged vertically upward into a uniform quiescent ambient. 
Temperature profiles were measured by a rake of thermistors in the flow 
field. Of interest in these experiments was the presentation of data 
which showed temperature residuals as measured near the terminal point 
of rise at the surface. Froude numbers of 1 to 50 and nozzle diameters 
of 0.425 and 0.55 inches were used. 

A summary of the range of parameters studied in various experiments 
is presented in Table (T3-4). 
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Table (T3-4). Summary of Experimental Parameters 
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IV. EQUATIONS OF STATE AND AMBIENT STRATIFICATION MODELLING 



The determination of the discharge Froude number and the buoyancy 
force requires the evaluation of the temperature, t, and concentration, 
c, effects on density. Often the density, p(t,c,p) is a sufficiently 
linear function of both t and c over the range of temperature and concen 
tration difference between the jet and the ambient. Then density differ 
ences may be accurately estimated in terms of the two volumetric 
coefficients of expansion. 



3 = 






and 



c,p 



_L(3£.) 

p r 3C t,p 



where p r is some reference value of p, say pg = p(tg,Cg,Pg). Then 



Ap = -Bp r At and Ap = -yp r Ac 



for the separate t and c effects on density. The density at some t and 
c, in terms of the initial jet density P g, is written as 



p(t,C,Pg) ~ p(tg,Cg,Pg) ” P ( t g , C g , P g ) [ 3 ( t “ t Q ) “ Y ( C ” C g ) ] 



or 



P = Pg [1 - 8(t-tg) - Y(C-Cg)] 

In particular, the initial density difference between the local 
ambient and the initial jet, which appears in the Froude number, is 
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p aO " p 0 



B(t 0 -t a ) + r(c 0 -c a ) 



p 0 

In an unstratified ambient medium, the reduced, non-dimensional forms 
of equations (3-3,4 and 5) will contain the following three non-dimensional 
terms : 



0 -Q 

p aO p 0 



t — t 

(J!! §_) 

v t -t '* 

^0 l aO 



c -c 

(— — — ) 

V c aO 



With the values of X in Equations (2-2 and 3) the same, profiles of 
t and c are the same and 



t -t c -c 

/ rn a \ _ ,_m a 

t-t c-c 

^0 r a0 c 0 c a0 



) 



at all points along a trajectory at which (tg-t^) f 0 and ( c g _c a g) t 0- 
Thus, for constant values of s and y in an unstratified ambient. 



(!i_jn_) = ( 
p a0" p 0 



t -t 
_m a_ 

t 0" t a0 



c -c 

) = (r-r~) 

C 0 c a0 



and defined values of 3 and y are not required for integration, computa- 
tion, and solution, beyond initial definition of discharge Froude number. 

A similar situation exists for some special cases of linearly strati- 
fied ambients. If p = p(t), and density stratification is defined solely 
by a temperature stratification parameter 



3t a 

<^)D 



l 3Z 

It 



aO'V 
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Then the definition of 6 beyond the initial computation of discharge 
Froude number is not required. For such flows, 



a~ p m \ 
p a0" p 0 



t -t 
■ m 

v* 



a 

aO 



) 



pertains throughout. In addition, if p f p(c), an assumption that might 
be made in a jet with a tracer dye, the concentration and vertical momen- 
tum equations are not coupled. Concentration computations can proceed 
independently, even to the point of specifying a concentration stratifi- 
cation if desired. 

Parallel reasoning holds true for the case of p = p(c), in which case 
a density dependence on concentration stratification parameter would 
couple the concentration and vertical momentum equation. Presumably 
temperature would be uniform throughout the system in such a case. 

These simple, degenerate cases of the overall modelling problem are 
important because they represent conditions under which measurements 
are often taken. Modellers then specify these conditions to compare 
data with their analytical models. 

The underlying assumptions and limitations of such formulation are: 

(1) With no equation. of state incorporated in the computational 
process, only unstratified ambients may be acommodated if 

p = p(t,c), and then only if s and y are assumed constant. 

(2) If no equation of state is included in the model, density 
stratification in the ambient may be specified as either 



at 

^aO'V 



p = p(t) 
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or 




or 




t- — y , (directly). 

^ p aO~°o' 



Ambients with more than one density-constituent gradient cannot be 
accommodated, since temperature, concentration, and density residuals 
would be independent of each other and profiles could be dissimilar. 

The more general case of an ambient medium with temperature and 
concentration stratification requires an equation of state for solution 
of the governing equations. Such a temperature/ concentration/density 
relation may be used in one of two ways: 

(1) Internally in the calculational scheme, with ambient temperature 
and concentration gradients specified. The local density differ- 
ences are computed downstream from jet temperature and concentration 
decays calculated from the energy and concentration equations. 

(2) External to the actual integration calculations, by using speci- 
fied ambient stratification to calculate the ambient density 
variation a priori. The resulting ambient density gradient is then 
used in an equation of density excess or deficiency of the form: 



HS^O Jq u(p a" p m )rdrd6} 



do 2tt °° 



35 ^/ / Urdrd<j>} 



0 0 
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assuming the same Gaussian profile for density as for two of 
its constituent properties, temperature and concentration. 

This extra equation effectively uncouples the temperature and 
concentration equations from the vertical momentum equation. 

Method (1) above was used by Hirst [15] . The second method above was 
used by Fan [9] . Other models specify the method of use of an equation 
of state, only as noted below. 

Hirst specified values of 3 and y, assumed constant throughout the 
downstream integration of the governing equations. The equation of 
state was: 




+ Y(c-C a0 ) 



In reality, of course, density of water is a complex and, under some 
conditions, highly non-linear, function of three variables--temperature, 
species concentration, and pressure. Many modelling situations might 
arise in which the range of these variables clearly argue against the use 
of constant values of y and 3, as well as the omission of the pressure 
dependence of density. Table (T4-1 ) shows some of the variation of 3 
encountered over small temperature ranges, for both saline water at 35 °h Q 
(typical of sea water) and fresh water. The pressure in both cases is 
one atmosphere. 

In their model, Riester, Bajura, and Schwartz [23] attempted to 
account for the effect of variable s by representing it in the form 



= a-| + a 2 t + a^t 



(4-1) 
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Table T4-1. Variation of 6 over limited temperature range, 
for pure and saline water 
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where the coefficients a-j , a£, an d a 3 were evaluated for the temperature 
range pertinent to a given calculation. The discharge Froude number was 
then "adjusted" using the relation 



<9°S a (V‘a» 1/2 



They summarized the need for an effective equation of state as follows: 

...if the initial temperature difference between the jet and 
the ambient is greater than 3°C, the effect of the temperature- 
density relationship of water on the trajectory must be con- 
sidered. This result is of particular significance when 
laboratory models are operated at high temperature differences 
to simulate prototype flow characterized by low temperature 
differences. .. .the effects of temperature on the trajectory 
can be minimized by using the (adjusted) Froude number ... 
to define flow conditions ... .Additional work is required to 
determine the correlation application to jet flows which 
experience reversing buoyancy for cases near 4°C. 

Salt water jets cannot simulate the characteristics of 
fresh water jets operating over large temperature differences 
since the salt models do not duplicate the temperature- 
density characteristics of the fresh water system. [24] 

In view of this last observation, it is noted that, although the proposed 



model attempted to compensate for the variable temperature effects on 
density, no dependence on salinity was included except for that which 
may have been implicit in the choice of the three coefficients in 
Equation (4-1). 

The ideal equation of state for buoyant jet modelling would compre- 
hensively represent density as a function of temperature, species concen- 
tration, and pressure. Such an equation should be applicable in all 
temperature, concentration, and pressure ranges of interest. Also, for 
most modelling of fresh/salt water systems, the concentration variable 
should be salinity. An equation of state which seems to fulfill these 
requirements will be suggested in the following section. 
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V. PRESENT METHOD 



The foregoing modelling system assumptions regarding the flow field 
also apply to all the calculations done in this work. Generally, they 
assume a steady, axisymmetric jet with negligible molecular transport, 
negligible streamwise turbulent transport, and small curvature effects, 
operating in a hydrostatic pressure field under the Boussinesq approximation. 
The governing equations in the zone of established flow follow. 

Continuity of mass equates the downstream change in total mass of 
the jet to the mass of fluid entrained. The entrainment concept results 
in the following: 

i{/ / pUrdrdcfi} = p (2iraU B) 

ai 0 0 m 

The first Boussinesq approximation, concerning density level, yields: 

2 °° 

7 U/ 77 / Urdrd<j>} = 2™U m B (5-1) 

ai o o m 

Horizontal momentum is conserved in a hydrostatic pressure field. 

Thus, the change in horizontal momentum within the jet is equal to the 
horizontal momentum of the fluid entrained: 

to I 27 " / P U 2 coserdrd<j,} = P UE 
ai o o 3 v 

Again, invoking the Boussinesq approximation: 
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( 5 - 2 ) 



cos erdrd<j>} 

0 0 



2tt 00 




The change in the vertical momentum of the jet system is the result 
of the action of the buoyancy force over the extent of the jet: 



In a stratified ambient, the ambient temperature, as well as the 
temperature of the jet, may be variable along the trajectory. The energy 
equation relates the change in energy of the system, as expressed by a 
temperature excess or deficiency relative to the ambient, to the rate 
at which the ambient temperature is changing. 



Assuming constant specific heat and the Boussinesq approximation. 



Similarly, the change in concentration excess relative to the 
ambient is related to the rate at which the ambient reference concentra- 
tion is changing: 



4^{/ / p sinerdr} = / / (p -p)grdrd<j> ( 5 - 3 ) 

0 0 0 0 




ii/ C P U(t-t )rdrd<|>} = - ^ / / pC Urdrd<t> 

05 0 p a 0 0 p 




Jk/ / U(t-t Jrdrd*} = -ST-I / Urdrdd, (5-4) 
0 0 3 ab 0 0 





Again, invoking the Boussinesq approximation: 
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oo 



2tt 00 dc 2tt 

/ U(c-c,)rdrd<j>}' = - -ts— j / Urdrd<j> (5-5) 

0 0 a 05 0 0 

Geometrically, the incremental trajectory of the jet may be 
described by: 



dX = 


dS cos 9 , 


(5-6) 


dZ = 


dS sin 9 


(5-7) 



It will be assumed that the Gaussian velocity, temperature, and 
concentration profiles cited previously are valid for the zone of flow 
establ ishment. 

Appendix A contains the development of the governing equations from 
their stated form above to their reduced form. In Appendix B the equa- 
tions are non-dimensional ized. The dimensionless equations appear in 
Table (T5-1). The forms of the equations for degenerate cases are 
presented as follows: 

(1) Flowing unstratified ambient--Table (T5-2); 

(2) Quiescent, stratified ambient--Table (T5-3); 

(3) Quiescent, unstratified ambient--Table (T5-4). 

The above equations apply only in the zone of established flow, 
which begins at S g , see Figures F-2-1 and F-2-2. Estimates of the end 
of the region of flow establishment have been made. The procedures 
followed here in calculations are based on the experimental work of 
Abraham [2] and the analytical development by Hirst [13]. They are also 
consistent with assumptions made in references [14], [15], and [12]. 

The initial conditions are: 
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Table (T5-1). Non-dimensional Differential Form of the Governing 
Equations for a Flowing, Stratified Ambient 
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Table (T5-2). Non-dimensional Differential Form of the Governing 
Equations for a Flowing, Unstratified Ambient 
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Table (T5-3). Non-dimensional Differential Form of the Governing 
Equations for a Quiescent, Stratified Ambient 
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VERTICAL TRAJECTORY ^( z ) = sin 



Table (T5-4). Non-dimensional Differential Form of the Governing 
Equations for a Quiescent, Unstratified Ambient 




47 



U = IL 
me 0 



X . = X = 1.16 

t c 



B e = D //2 



’V'.le (x t + '> 



nirVo 2xl 



<V‘a>e < x c +1 > 

TcITcTJ 



m a'O 



2x 



The value of S g , relative to the point of discharge, is 



S e /D = 6.2 F > 40 

S e /D = 3.9 + 0.057F 2 5 <_ F < 40 

S e /D = 2.075 + 0.425F 2 1 < F < 5 

S e /D =0. 0 < F < 1 

All that remains to be specified to integrate the equations down- 
stream along the S coordinate are the entrainment function to be used, 
the ambient stratification condition, and the density equation of state. 

The calculations used an IBM 360 and IBM 3033 computing system. 

The method of integration is a trapezoidal rule expression of the 
governing equations with adjustable step size. Two versions of the 
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general method were used here. The first is used to evaluate flows in 
unstratified ambients, or in ambients characterized by linear gradients 
in temperature, concentration, or density itself. No explicit equation 
of state is necessary, only 6 and y. The second version, with the same 
computational technique, uses an accurate density equation. 

Since the objective is to develop a predictive model applicable to 
large scale discharges in ocean water, the density equation used in 
version two is that of Gebhart and Mollendorf [11]. This relation 
correlates the temperature dependence of density as an expansion around 
the density extremum temperature at any given level of pressure and 
salinity. The resulting expression is fitted to a comprehensive set 
of experimental data over temperature, salinity and pressure maximums 
of t = 20°C, s = 40%, and p = 1000 bars absolute. Density is in kg/rn^. 
This equation is given in Appendix C. 

For each of these calculations, the density equation was first used 
to establish the density field in the ambient. This requires input of 
at least one reference temperature and salinity level and appropriate 
gradients, or specification of the complete temperature and salinity 
fields. The water surface pressure is assumed to be one atmosphere. The 
relation is first used at the surface, and density is calculated pro- 
gressively downward, using hydrostatically integrated values of pressure 
along the way. 

Once the ambient density field is established, the relation is then 
used to calculate the density terms in the vertical momentum equation at 
each calculational interval, using the jet temperature and concentration 
values predicted by their respective equations. Ambient density and 
pressure are assumed to be linear between the points calculated for the 
ambient. 
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The relation is also used to calculate the initial density differ- 
ence between the jet and the ambient at the point of discharge, given 
temperature and salinity level, or to calculate the initial temperature 
of the jet, given Froude number and salinity. 

The advantages of using an accurate density relation include: 

(1) The full density dependence on temperature, salinity and pressure 
is taken into account in all calculations. 

(2) The actual ambient density field of either a hypothetical or 
real circumstance is accurately established, utilizing appro- 
priate temperature and salinity input. 

(3) The relation is readily applicable to computerized analysis, and 
is computationally compact. 
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VI. RESULTS 



Comparison between the jet trajectories and properties calculated 
here and earlier calculations, using the same entrainment function, has 
been generally good. This is in spite of differing computing techniques 
used in different studies. The resulting trajectories will be compared 
with earlier calculations. 

A. UNSTRATIFIED, QUIESCENT AMBIENT 

This jet/ambient system represents by far the most frequently modelled 
case. Comparison here is for horizontally discharged buoyant jets. 

Figures (6-1) and (6-2) compare the trajectories predicted by these 
calculations, for F = 10 and 20, with the calculations of Abraham [3], 
both using a = 0.085. Correspondence between the two is consistent, with 
the present method predicting slightly less horizontal penetration of 
the ambient than Abraham's model at all Froude numbers compared. 

Figures (6-3) and (6-4) show the uniformly close correspondence between 
these calculations and those of Hirst [14], [15], both using Hirst's en- 
trainment function, entry 3 in Table T3-2, for F = 4,6,8 and 10. 

Figures (6-5) through (6-8) again compare present calculations with 
those of Hirst, but with both using the entrainment functions a = 0.057 
and a = 0.082. Data of Fan is included for comparison where available. 

The results, combined with those in Figures (6-3) and (6-4) indicate that 
within the Froude number range compared, correspondence of present calcu- 
lations with Hirst's model is very good regardless of the choice of 
entrainment function. 
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igure 6-1. Comparison of present method with model of Abraham [3] 
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Fioure 6-2. Comparison of present method with model of Abraham [3] 
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Figure 6-3. Comparison of present method with model of Hirst [13] 
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Figure 6-4. Comparison of present method with model of Hirst 
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Figure 6-5. Comparison of present method and model of Hirst, 
both using a = 0.057, and data of Fan [7] 
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Figure 6-6. Comparison of present method and model of Hirst, 
both using a = 0.057 
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Figure 6-7. Comparison of present method and model of Hirst, 
both using a = 0.082, and data of Fan 
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Figure 6-8. Comparison of present method and model of Hirst, 
both using a = 0.082 




Comparison of calculations with those of Fan [9], using a = 0.082, 
is shown in Figure (6-9). In this case, the present method predicts 
greater horizontal penetration of the jet and, therefore, somewhat less 
vertical penetration at most values of F. 

A comparison of the present calculations, using no equation of state, 

with the calculations of Riester et al . [23] , both using 

2 2 1/2 

ct = [(0.057 cose) + (0.082 sine) ] ' , showed reasonably good agreement 
for high temperature ambients, where the 3 correction factor of Riester 
et al . resulted in little adjustment of the Froude number. Progressively 
poorer comparison was evident with their prediction for lower temperature 
ambients, where the Froude number adjustment became larger. This is shown 
in Figure (6-10) . 

Figures (6-11) and (6-12) compare the calculated jet trajectories at 
F = 1,2,4,6,8,10,50,100,150, and 200 for the five quiescent ambient 
entrainment functions. The ambient is unstratified and no equation of 
state is used. The models which relate to plume and non-buoyant jet 
flows understandably begin to deviate significantly from the norm of the 
trajectories when F is not in their intended range of use. This is 
particularly true for the buoyant plume entrainment function, a = 0.082 
(model 2). It consistently predicts lower entrainment and higher trajectory 
at all but the lowest Froude numbers. Albertson's momentum jet function, 
a = 0.057, (model 1) consistently predicts the highest entrainment and, 
therefore flattest, trajectory at all values of F. The "end-point" 
correlation of Riester et al . , model 5, is also seen to predict higher 
entrainment than the norm for all values of F. 

The generally accepted values of entrainment functions for the extrema 
of Froude numbers are a = 0.082 for F = 0, and a = 0.057 for F = ». The 
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Figure 6-9. Comparison of present method with model of Fan [9] 
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Figure 6-10. Comparison of present method with model and 
data of Riester et al . [23] 
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Figure 6-11. Comparison of entrainment models, quiescent ambient medium 
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Figure 6-12. Comparison of entrainment models, quiescent ambient medium 



two models which most closely reflect the Froude number/entrainment 
relation over the range of F considered here are the function of Hirst 
(model 3) and the empirically derived function of Davis et al . (model 4). 
Hirst's function actually predicts the most buoyant, least entraining 
behavior for the lowest Froude numbers (<_ 10). 

All models predict strikingly similar behavior in the decay of 
centerline velocity, temperature, and concentration. The decay of jet 
properties along the pathlength S are not strongly dependent on F or 
entrainment function. Figure (6-13) shows the superposition of center- 
line velocity decay, U/Uq, for all jets depicted in Figure (6-11) (F = 10, 
50,100,150,200). The trend is seen to be an extremely rapid initial 
velocity decrease immediately after discharge. For jets with F >_ 50, 

U/Uq =0.1 around 50 diameters along the trajectory. Further downstream, 
the residual velocities decrease less rapidly. Figure (6-13) also shows 
that higher velocities persist downstream at smaller values of F, i.e., 
for more buoyant or less vigorous jets. 

Close inspection of Figure (6-13) reveals an apparent anomaly. For 
F = 10, the Hirst model predicts the lowest residual velocity along most 
of the trajectory. The highest residual velocity is predicted for the 
buoyant plume entrainment function a = 0.082. Yet Figure (6-11) for 
F = 10 shows that these two entrainment functions predict the highest, 
most buoyant trajectory, with the Hirst function representing the extreme. 
The behavior of the Hirst function, on one hand displaying buoyant charac- 
teristics indicative of low levels of entrainment, and on the other 
hand displaying high velocity decay indicative of high entrainment, 
seems to contradict the consistent behavior of the buoyant plume function. 



65 



LO 




o 

o 

CM 

o 

4-> 

o 



oo 

4-> UO 
OJ 

•r-3 O 
4-> 
UO 

CM r— 



S- 

o 



>> 

03 

U 

CD 

*o 



o 

o 



<D 



00 






( V 
S- 
13 



•n/n 



66 



entrainment models 




This behavior may be explained by close examination of Figure (6-14), 
an enlargement of the F = 10 results on (6-11). In the first 20 diameters 
of trajectory, the Hirst model is seen to be lower because of more rapid 
entrainment in comparison to the buoyant plume model. It then curves 
rapidly upward, displaying the most buoyant trajectory farther downstream. 
The velocity level apparently decreased early due to rapid entrainment. 
However, the jet still retained enough buoyancy to cause the sharp upturn 
characteristic of the later trajectory. These results are indicative of 
the importance of the early level of entrainment in determining properties 
much further downstream. 

An even more uniform behavior among the group of 25 calculations con- 
sidered in the illustration of velocity decay is shown in a superposition 
of temperature difference decay curves in Figure (6-15). Jets in this 
F range of 10 to 200 have temperature residuals of less than 0.2 at 20 
diameters downstream and less than 0.1 at 50 diameters. Both of these 
decay curves are in close agreement with existing data. Note that in 
this formulation, concentration residuals decay in the same way. 

An example of the radial or diameter growth of a jet, using entrain- 
ment model 2, is shown in Figures (6-16) and (6-17). These figures 
illustrate that the physical extent of a jet generally increases more 
rapidly with increasing F. 

In summary, the various past entrainment models for jets in a quies- 
cent ambient medium predict reasonably similar behavior. Most notable 
exceptions to this conclusion are the results using the buoyant plume 
entrainment function, model 2, at higher F. This is well outside the 
range of intended use. Further, all of the models are in reasonably 
close agreement with meager existing data for small diameter jets. 
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Figure 6-14. Comparison of entrainment models, quiescent ambient medi 
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= 10 to 200, entrainment models 
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Figure 6-16. Radial growth of buoyant momentum jets, quiescent ambient medium 
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Figure 6-17. Radial growth of buoyant momentum jets, quiescent ambient medium 




Twenty years of "entrainment function generation" in the literature not 
withstanding, it may be stated that in a quiescent ambient, predicted 
jet behavior is not strongly dependent on the particular entrainment 
model used. 

B. FLOWING UNSTRATIFIED AMBIENT 

The two entrainment functions proposed for a buoyant momentum jet in 
a flowing ambient are models 8 and 9. It will be seen that they do not 
result in agreement between predicted results characteristic of the 
quiescent ambient models. 

Figure (6-18) shows the trajectory predicted by the two entrainment 
functions for a vertical jet in a horizontal cross flow. The conditions 
are R = 0.25 and F = 2.83, in an unstratified ambient. This set of values 
of R and F is especially significant because it is one of the few cases 
in which the models are known to be in close agreement. For all F less 
than 2.83, model 9 predicts a higher, more buoyant trajectory than model 
8. Conversely, at F > 2.83, model 9 generally predicts lower trajectory 
and somewhat more rapid decay of velocity, temperature, and concentration. 
For F = 2.83, model 9 predicts lower trajectory for all R < 0.25, and vice 
versa. Figure (6-19) shows the effect of lessening the ambient flow 
rate to R = 0.125 for a jet of F = 2.83 with all other conditions the 
same as Figure (6-18) . 

Figures (6-20), (6-21) and (6-22) illustrate that the disparity between 
the trajectory predictions of the flowing ambient entrainment models is 
not limited to the case of discharge nromal to a cross flow. In these 
cases, the discharge angle is 45°. The same behavior is evident--reason- 
ably good agreement for F = 2.83 and R = 0.25 as in Fig. 6-20. However 
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Figure 6-18. Comparison of flowing ambient entrainment models 
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Figure 6-19. Comparison of flowing ambient entrainment models 
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Figure 6-20. Comparison of flowing ambient entrainment models 
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Figure 6-21. Comparison of flowing ambient entrainment models 
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Figure 6-22. Comparison of flowing ambient entrainment models 



model 8 predicts a higher trajectory for R > 0.25, Fig. (6-21), and 
model 9 a higher trajectory for R < 0.25, Fig. (6-22). 

As F and R are increased, the upward penetration of the jet decreases. 

This is due to more rapid entrainment and lower initial buoyancy of the 
high Froude number jets, as well as the stronger horizontal "sweeping" 
effect of the more rapidly flowing ambients. In most cases, however, 
the increase in these effects do not lessen the disparity between the 
relative magnitudes of rise predicted by the two entrainment models. 

Figure (6-23) depicts another case of relatively good agreement at F = 10, 

R = 0.125, 6q = 90°, compared with the data of Fan at these values. How- 
ever, Figure (6-24) shows that model 8 diverges from the other two indi- 
cated trajectories when R is increased to 0.25 at the same value of F. 

Figure (6-25) shows the same disagreement between models 8 and 9 at 0g = 45°. 
For higher Froude numbers, the disparity persists, as shown in Figures 
(6-26) to (6-31). 

The results for initial pure "coflow", in which the jet is discharged 
horizontally parallel to the ambient flow, 9 q = 0, are shown in Figures 
(6-32) to (6-37). The conditions are F = 20 and 40 for R = 0.125, 0.25, 

0.4 and 0.8. The same characteristic disagreement is evident. At small 
values of R, the disagreement between the models is 15 to 40%. At larger 
R values, the disagreement approaches 100%. The two models are inconsistent 
throughout. These high R values are beyond the data used in the fit by 
Ginsberg and Ades in model 9. 

In summary, the disparity between the predicted jet trajectory using 
the two flowing ambient entrainment functions is generally most apparent 
in flows with significant vertical components in their trajectory. Obvi- 
ously, the models are very sensitive to how the cross flow entrainment 
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Figure 6-23. Comparison of flowing ambient entrainment 
models and the data of Fan 
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Figure 6-24. Comparison of flowing ambient entrainment models 
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Figure 6-25. Comparison of flowing ambient entrainment models 
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Figure 6-26. Comparison of flowing ambient entrainment models 



o 




a/2 



83 



Figure 6-27. Comparison of flowing ambient entrainment models 
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Figure 6-28. Comparison of flowing ambient entrainment models 
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Figure 6-29. Comparison of flowing ambient entrainment models 
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Figure 6-30. Comparison of flowing ambient entrainment models 
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Figure 6-31. Comparison of flowing ambient entrainment models 
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Figure 6-32. Comparison of flowing ambient entrainment 
models for case of pure coflow 
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Figure 6-33. Comparison of flowing ambient entrainment 
models for case of pure coflow 
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Figure 6-34. Comparison of flowing ambient entrainment 
models for case of pure coflow 
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Figure 6-35. Comparison of flowing ambient entrainment 
models for case of pure coflow 
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Figure 6 - 36 . Comparison of flowing ambient entrainment 
models for case of pure coflow 
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Figure 6-37. Comparison of flowing ambient entrainment 
models for case of pure coflow 



effect is accounted for. These low Froude number, low R flows often 
result in calculated differences of vertical rise of many diameters in 
magnitude. This disagreement is found at all discharge angles. At 
higher values of F and R, the absolute difference in predicted vertical 
rise becomes small, but relative to the total rise of the jet, the dif- 
ference may be many times greater than for lower values of F and R. 

C. STRATIFIED AMBIENTS 

A stable stratified ambient, in which density increases with increasing 
depth, has the general effect of restricting the vertical motion of an 
initially buoyant jet. This is because of the compound effect of strati- 
fication and of the jet losing buoyancy due to entrainment of ambient 
fluid. This occurs regardless of the direction of buoyancy, ambient 
stratification, or of the magnitude of the stratification. 

Figure (6-38) illustrates the effect of increasingly strong stable 
stratification, due to an ambient temperature gradient. The calculation 
is for an upwardly buoyant jet of F = 50, discharged horizontally into a 
quiescent ambient. For the least stratification. A, the jet has followed 
a trajectory similar to the unstratified case, replotted from Figure 6-11. 
There is a very slight vertical restriction of the trajectory. In case 

B, near the end of the computed trajectory, a noticeable change in trend, 

or re-bending of the trajectory, is evident. The point of inflection in 
this curve is very significant. Calculations show that this is the verti- 
cal level at which the jet has become neutrally buoyant, due to entrainment 
on one hand and the decreasing density of the surrounding fluid on the 

other. Further upward rise beyond this point is due solely to the momentum 

then existing in the jet. This momentum is gradually decreased by the 
downward force of negative buoyancy. 
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Figure 6-38. Buoyant jet behavior in a stratified ambient 



Case (C) shows the eventual course of events for a jet which has 
become vertically "trapped", due to increased ambient stratification. 

The trajectory has undergone complete re-bending, to horizontal flow 
with no buoyancy. Standard practice in entrainment modelling, which 
will be followed here, is to halt computation at the point of maximum 
rise. Beyond this point, all but the most vigorous jets will have 
acquired "far field" characteristics. Assumptions of axisymmetric jet 
shape, tenuous at best in a stratified ambient, certainly cannot be made 
beyond this point, and initial momentum and buoyancy have become largely 
irrelevant to future behavior in most cases. 

Figures (6-39), (6-40), and (6-41) illustrate the effect of the 
same degrees of stratification on jets at a higher discharge Froude number, 
F = 100, 150 and 200. The general behavior characteristics are similar 
to those previously noted, with ambient stratification further restricting 
the vertical rise already limited by low initial buoyancy. 

The intuitive reasoning that the trajectory of initially more buoyant 
jets is less affected by stratification is supported by these calculations. 
Figures (6-42), (6-43), and (6-44) compare the trajectories of various 
Froude number discharges in an environment of three different stratifica- 
tion levels. The lowest F jets are seen to be negligibly affected over 
this range of stratification, while the higher F jets exhibit varying 
degrees of trajectory flattening, vertical trapping, and terminal rise. 
Given strong enough stratification in the ambient, of course, any initially 
buoyant jet will eventually experience negative buoyancy and terminate its 
rise at some elevation in a sufficiently extensive ambient. No elevation 
overshoot is seen. 
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Figure 6-39. Jet behavior in a stratified ambient 
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Figure 6-40. Jet behavior in a stratified ambient 
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Figure 6-41. Jet behavior in a stratified ambient 
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Figure 6-42. Effect of stratification on jet behavior 
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Figure 6-43. Effect of stratification on jet behavior 
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Figure 6-44. Effect of stratification on jet behavior 



The above circumstances demonstrated are typical of past jet modelling 
in a stratified ambient. Stratification, as indicated in the earlier 
discussion of equations of state, is usually represented in terms of a 
single, constant temperature, salinity, or density gradient. Rarely, 
as in Reference [15], temperature and salinity gradients may be treated 
together. 

The capability of a comprehensive equation of state to more realis- 
tically describe a particular water ambient and the behavior of a jet 
within it is illustrated in Figures (6-45) to (6=-50). Figure (6-45) 
characterizes the temperature, salinity, and density stratification of 
an area of the Northern Pacific Ocean to a depth of 500 m. June tempera- 
ture and salinity data for a point 53° 04' N, 175° 35' W were taken from 
Reference [4], and approximated by a series of 5 temperature and 4 salinity 
gradients. The Gebhart-Mollendorf relation was then used to establish 
the corresponding density field shown, which includes the contribution 
of hydrostatic pressure. 

Figure (6-46) shows the predicted trajectories of 5 buoyant momentum 
jets, for F = 10,20,30,40 and 50, discharged horizontally in this density 
field at a depth of 200 m. In each case, the initial jet diameter is 1 m. , 
while F ranges from 10 to 50. Initial salinity difference between the 
jet and ambient was assumed to be zero. Buoyancy was due to elevated 
temperature, however the effect of the salinity stratification gradient 
is included in the calculation of jet density and buoyancy along the 
trajectory. 

All five of these jets exhibit the progression from positive to nega- 
tive buoyancy, as well as a level of maximum upward penetration. All dis- 
play a rather smooth curvature and re-curvature, due to the relatively 
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Figure 6-45. Temperature, salinity and density for Northern 
Pacific Ocean, 53° 04' N, 175° 35' W, June 
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Figure 6-46. Trajectory of 5 jets discharged into the ambient of Figure 6-45 
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Figure 6-47. Temperature, salinity, and density profiles for tropical 
Atlantic Ocean, 02° 03' S, 39 °20' W, February. 
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Figure 6-48. Trajectory of 5 jets discharged into the ambient of figure 6-47 
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Ficjure 6-49. Temperature, salinity, and density profiles for 
Arctic Ocean, 81° 28' N, 8° 05' E, October 
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Figure 6-50. Trajectories of 5 jets discharged into the ambient of figure 6- 49 



constant density gradient encountered in the vertical range of penetra- 
tion. Again, no overshoot occurs. 

Figure (6-47) illustrates another characteristic oceanic temperature- 
salinity-density field, in this case from the tropical Atlantic at a point 
02° 03' S, 39° 20' W. The data for this construction was taken in Febru- 
ary and is from Reference [21]. The trajectory of jets F = 10,20,30,40 
and 50 discharged into this density field at a depth of 200 m. is shown 
in Figure (6-48). The vertical penetration of these jets in this case 
is significantly less than calculated for the Northern Pacific density 
profile, due to the much stronger density gradient in the tropical ocean. 

An October temperature-sal inity-density field for the Arctic Ocean 
at 81° 28' N, 8° 05' E is shown in Figure (6-49). Data is from Reference 
[17]. The strong salinity discontinuity near the surface is due to melt 
from the adjacent Arctic ice pack. The behavior of jets discharged into 
this ambient, Figure (6-53), is especially interesting compared to Figures 
(6-49) and (6-51), because of the varying recurvature rates encountered. 

For F = 30, the jets exhibit a smooth curvature and recurvature. The 
F = 10 and F = 20 jets, however, show a much sharper recurvature because 
of the encounter with the steeper density gradient above the 100 m. level. 

These examples illustrate the advantage of using a full equation of 
state in entrainment modelling. The dependence of the model on assumed 
values of s and y is eliminated, as is the requirement to express ambient 
stratification in terms of constant temperature, salinity, or density 
gradients. The accuracy of the ambient density characterization is limited 
only by the availability and spacing of the field data. Ultimately, the 
use of a full equation of state can reduce one of the primary uncertain- 
ties in entrainment model 1 ing--that of evaluating the temperature- 
salinity-pressure-density relationship. 

no 



VII. EXTRAPOLATION OF PRESENT METHODS TO PARAMETERS OF INTEREST 



In reviewing the progress to date in entrainment modelling of buoy- 
ant momentum jets, some rather striking matters become apparent. The 
first of these is the very small physical scale of the experimental studies 
which has been made to underlie subsequent modelling schemes. Jet dis- 
charge diameters were mostly about 1 cm, the largest being 3.9 cm. The entire 
jet trajectories were usually only a few meters in length. Great lengths 
arose in some of the experiments done with flowing ambients, but even 
these studies used similarly small diameter discharges. These limitations, 
of course, arose, at least in part, through the innate limitation in size 
and economy which applies to experimental work. 

Still, the data base from these small scale discharges are often used 
to develop models for real world jets, such as power plant discharges and 
sewer outfalls. These are commonly hundreds to thousands of times larger. 

Such upscaling of the implications of the experiments implies a faith in 
the similarity properties of the densimetric Froude number. This may or 
may not be justified. Further, such a large scale difference between 
supporting data and real system projection, in the case of an underwater 
momentum jet, does not necessarily imply similar scaling in turbulence 
within the jet and in the immediately surrounding ambient. Fan [7] found 
that the effect of ambient turbulence on jet behavior was profound. It 
is not unreasonable to expect that the scale of turbulence within the 
jet, or the relative scale between ambient turbulence and jet dimensions, 
may not have equally important effects. 

The second striking aspect of experimental and analytical work is the 
relatively low densimetric Froude number range used. Admittedly, more 
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buoyant, lower Froude number jets are more "interesting" in terms of 
possible trajectories, and may be of more importance in many environmen- 
tal modelling situations. Froude numbers in the range of 10 to 40 are 
common in the literature. Values as high as 100 are rare. There are a 
number of "real" cases, however, in which entrainment modelling might be 
used, where Froude numbers might be in this higher range. Very basic 
calculations dealing with unclassified physical and propulsive properties 
of submarines, for instance, indicate that the propeller outflow, with 
all condenser efflux entrained and mixed within, might represent a Froude 
number in the range of 50 to 200. Ship or submarine condenser efflux by 
itself might be expected to have a range inclusive of all Froude numbers 
to about 200, depending on discharge size, flow rate, and thermal loading. 

A third limitation in the present state of entrainment modelling is 
the inconsistent predictive quality of the collection of available flowing 
ambient entrainment functions in their supposed range of applicability. 

A fourth shortcoming is the range of R in which data is available to guide 
modelling. Very little experimental data exists for R greater than 0.25. 
Almost no data exist for an initially co-flowing ambient. The data is not 
sufficient to definitively select either a general entrainment function 
or even one valid over a restricted range of F and R. 

Given the limitations of the current background data to support any 
modelling technique, a high level of confidence in predicting large scale 
jet trajectories, involving conditions of ambient flow, turbulence, or 
high Froude number discharges is unwarranted. There are simply too many 
unknowns to proceed confidently, limited by contemporary methods and 
information. 

Toward the objective of developing large scale entrainment modelling, 
the following suggestions are made. 
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VIII. RECOMMENDATIONS 



If entrainment modelling is to be made applicable to large scale 
discharges involving relatively high ambient flow velocities, a more 
comprehensive data base than now exists needs to be established. Such 
a data base would include: 

(1) higher Froude number flows, up to a range of 200. 

(2) larger scale discharges, with particular attention to turbulence 
scale and its effects. 

(3) higher R coflows. 

Future analytical models, drawing from such an enhanced data base, 
need to specifically address: 

(1) determination of a dependable entrainment function for submerged 
jets discharged to a flowing ambient. 

(2) determination of methods to include turbulence effects in the 
model . 

(3) inclusion of an accurate and comprehensive equation of state, 
such as the one proposed here, in the model. 
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APPENDIX A 



A. ASSUMED GAUSSIAN DISTRIBUTIONS 



U - U m exp[.^-] 



T-T _ 2 

T^T~ = exp[ 72l2 ] 

m a X B 



C-C. 



C -C 
m a 



exp [ p p ] 

X 



p a~ p 

p a‘ D m 



exp [-^r] 

X 



B. DEVELOPMENT OF THE GOVERNING EQUATIONS 
1 . Continuity 



dS^o f Q Urdrd * } = 2 ™ u m B 



^2tt/ Urdr} = 27TaU m B 



Substituting Equation (2-1) 



3T : L U m expl^rdr) - aiyi 

U D 



( 2 - 1 ) 



( 2 - 2 ) 



(2-3) 



(5-1) 
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d B 2 

^-{U V = all B 
db m 2 m 



4^{ll B 2 } = 2aU B 

dS m m 



(A-l ) 



For the case of a coflowing ambient, the entrainment rate is 
expressed as a function of jet centerline velocity and jet width which 
is more complex than the simple linear relationship of the quiescent 
ambient case. For coflowing ambients, 



2tt » 

TnKj / Urdrd*} = 2irf(a ,11 , B,R) 
° 0 0 



By development identical to the quiescent ambient case, 



3s< U ra B 2> * 2f <“’ U ,n’ B > R) 



(A- 2) 



2. Conservation of Energy 



2tt °° 

/ / U(t-t )rdrd<|>} 

db 0 0 



dt. 2 tt 

r 

0 0 



/ / Urdrd<}>} 



(5-4) 



Integrating and substituting Equations (2-1) and (2-2), 



d? 2 ” U m ex P^ ( V t a )exp[ - 27 lrd, ' ) 



dt oo 2 

1 f :2 " U ra exp(^-]rdr; 
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oo 



dS^ U m (t m" t a ) / Q exp[ } ]rdr} 



dt a /° -r 2 

“dS^m f Q ^P^ ]rdr} 



2 2 

4^u m (t m -tj[-^ — ]> = 

dS m m a 2 (a 2 +!) 



dt R 2 
dS^m 2 i 



(A-3) 



3. Conservation of Concentration or Other Scalar Species 



jci / / U( c— c ) rdrdcf> } 

0 0 



dc 2tt oo 

a [f / Urdrd* } 



dS 



(5-5) 



0 0 



By development identical to conservation of energy in (2) above, using 
Equation (2-3), 




U m ( V c a>'- 2 r~ li 

m m a 2 (a 2 +1) 



dc 

a 

dS 



R 2 

V 



(A-4) 



4. Conservation of Horizontal Momentum 
For the case of the quiescent ambient. 





co 




cos erdrdd)} 



0 



(5-2) 




0 



Substituting Equation (2-1), 
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0 



2ttU^cos 9 / exp[ ~ ^ ~ 3rdr} 
0 B 



^{27rU^cose[i ^-]} = 0 



-r^CU^B COS0} - 0 

db m 



For the case of the flowing ambient. 



H ‘■TT °° O 

/ / Ucos0rdrd<j>} = 2Tif(a,U m ,B,R)U a 

^2, j U 2 cos0rdr} = 27rf(a,U m ,B,R)U a 

Substituting Equation (2-1), 

°° 2 

^j{U 2 cos0 / exp[— j— ] rdr} = f(a,U m ,B,R)U a 
0 B 

2 

^ u ‘ c °se4V } = f(c,U m ,B,R>U a 

l^coseB 2 } = 4f(a,U in ,B,R)U a 

5. Conservation of Vertical Momentum 

. 2TT 00 r) 2iT °° 

/ olrsin0rdrd<|>} - J f (p -p)grdrd4> 
0 0 0 ' 0 a 



(A- 5) 



(4-6) 



(5-3) 
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00 00 p -p 

^-[2tt / U 2 sin0rdr} = 2n f (— — )grdr 



Substituting Equations (2-1) and (2-2), 



^sine / q exp[^f-]rdr} = ^ 9 / Q exp[^]rdr 



HS^ U m Sin0 ^2 B 2 ]} = p 



P a" P m g( i¥ ) 



^CU 2 sineB 2 } = ^-^(2gA 2 B) 
ab m p 



(A-7) 



6. Horizontal Trajectory 



dX = dX cos 0 



— X 
dS A 



= COS 0 



(A-8) 



7. Vertical Trajectory 



dZ = dS sin 0 



dS 1 “ sin 9 



(A-9) 
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APPENDIX B 



A. DIMENSIONLESS VARIABLES 

The following dimensionless quantities will be used in the non- 
dimensional ization of the governing equations: 



u = 



m 



m 



b = nr- 



s = 



X = 



z = 



= 




At o 


1 vw 


AC a 


(c a2- c al> 


t> 

O 

o 

1 


(c 0' c a0 ) 


At m 


<V‘a> 


At 0 


(V‘aO^ 


Ac m 


(V c a> 


At 0 


''o-'ao' 



B. NONDIMENS IONAL IZATION OF THE GOVERNING EQUATIONS 
1 . Conti nuity 

For the case of the quiescent ambient. 




u ,/> 



2ctU B 
m 



( A- 1 ) 



Substituting, 
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n-HU U n b 2 D 2 } = 2au U n bD 
as D m 0 mu 

iriu b 2 } = 2aU b (B-l ) 

as m m 

For the case of the flowing ambient, 

^nl) m B 2 } = 2f(a,U m ,B,R) (A-2) 



f(a,U m ,B,R) = (0.057 +^sine)B( |U m -U a cose| +a 3 U a sine) 

^Uu B 2 } = 2(0.057 +^-sine)B(|U -U cose I +a,U sine) 

abm 1 m a o a 

Substituting, 



2 . 



4 - ft -[ u U b 2 D 2 } 
ds D m m 



= 2(0.057 +-y 



97 



di^u^h j = 2(0.057+- 



Conservation of Energy 



sine)bD ( I u m U m’ RU O cos0 | + a 3 RU 0 sine) 
■^sine)b( |u m ~Rcose |+a 3 Rsine) (B-2) 






2( a 2 +1 ) 



dt a U B 

Lr_£! \ 

liT 1 2 * 



(A- 3) 



Substituting, 
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At, 



H U m^n a+ (t n -t n ) 



X 2 b 2 D 2 



ds D m 0 At g 0 dO 2 ( ' ^+T ) 



4t a ^ t O‘ t aO^ 



b 2 D 2 



Atg ds D 



{u U 



m 0 2 



h At ,2.2 

A { U in. A_J> > 

ds m4t O (X 2 +l) 



At 



a {u b 2 } 



Atg ds m 



(B-3) 



3. Conservation of Concentration or Other Scalar Species 



^7 } ■ 



dc_ R 2 
§ 4 - y 

ds " u m 2 J 



(A-4) 



By development identical to conservation of energy in (2) above, 



A c , 2 , 2 
m a b 



AC 



Hs^mZcT 7T? 



-} = - 



0 (a +1) 



3 (LI b 2 > 



A c g ds' m 



(B-4) 



4. Conservation of Horizontal Momentum 



For the case of the quiescent ambient. 



|KU 2 B 2 COS9> = 0 



(A- 5) 



Substituting, 



^u 2 l) 2 b 2 D 2 cose} = 0 



^-{uVcose} = 0 



(B-5) 
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For the case of the flowing ambient, 



^U^coseB 2 } = 4f(a,U m ,B,R)U a (A-6) 

f (ct , U^B , R) = (0.057+^sine)B( |U m -U a cose |+a 3 U a sine) 

^-{U 2 coseB 2 } = 

4U a (O.O57+^sin0)B( l u m -U a cose |+a 3 U a sine) 



Subsituting, 

4-D-' u m u 0 b2 ° 2cos6) * 

4RU o (O.O57+^-sin0)bD ( | u^Ug-RUgCOse | +a 3 RUgSi n0 ) 

^j{u 2 b 2 cos0} = 4R(O.O57+-y^sin0)b( |u |T| -Rcos0|+a 3 Rsin0) 

(B-6) 

5. Conservation of Vertical Momentum 

^U 2 sin0B 2 } = ^^(2gx 2 B) (A-7) 

The denominator of the density term on the right side of Equation 
(A-7) was transposed from the Bousinenesq term on the left side of the 
equation. It will be assigned a value equal to the reference discharge 
density. Substituting dimensionless terms, 
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d f ..2,,2 k 2 n 2 -4 -«, - p a~ p m /0 „, 2 u 2 n 2 x 

ds D lU m U O b D sine ^ p ( 2 ~ A b D ^ 



^u*t> 2 sin0} 



p a“ p m ( p aO' p mO ) 
p 0 p aCf p mO U 2 



o -p 9 2. 2 

am \ 2A b 



t;u 2 b 2 sine) = ( a ^ 

ds m PaO-°n.O F 2 



(B-7) 



6. Horizontal Trajectory 



dS X 



= cos e 



(A-8) 



Substituting, 



W~Q 



x D = cos e 



x = cos 0 



(B-8) 



7. Vertical Trajectory 






Z - sin 0 



(A-9 ) 



Substituti ng, 



ds D 



z D = sin 0 



ds 



z = sin 0 



(B-9 ) 
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